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Abstract 

We have recently calculated the second-order QCD corrections to the forward- 
backward asymmetry in e'^e~ annihilation. Here we recall the results and compare 
them to others in the literature. 



Experimental measurements of the forward-backward and left-right for- 
ward-backward asymmetries in e~^e~ annihilation to fermions provide some of the 
best determinations of the weak mixing angle sin^ Oejf ^ ■ In particular, the forward- 
backward asymmetry of b quarks is measured with a precision of about 2%, allowing 
an extraction of sin^ Oi^ff with almost per mille accuracy. However, since we are 
dealing with quarks in the final state, we must ensure that QCD corrections, both 
perturbative and non-perturbative, are understood to at least the same precision. 
For the perturbative corrections, this requires working to at least next-to-next-to- 
leading order (NNLO). 

To date there have been two 0{a'^) calculations, both in the massless ap- 
proximation and using a slightly different definition of the asymmetry than the 
experimental measurements, which use the thrust axis rather than the quark di- 
rection. The classic calculation of Altarelli and Lampe 0) determined the C(ag) 
coefficient numerically and found it to be small. This result has been the basis 
of all the experimental analyses since. However, the recent analytical calculation 
by Ravindran and van Neerven 0) obtained a coefficient about four times bigger. 
This discrepancy is comparable to the size of the experimental errors and needs to 
be resolved before the final electroweak fits to the LEPl data can be made. The 
C(ag)-calculation using the experimentally-used thrust axis definition, would also 
be highly desirable. 

We have recently performed a numerical calculation of the C(ag) correc- 
tions to the forward-backward asymmetry S). Anticipating the result, given below, 
we can say that to the precision required by experiment we confirm the result of 
Ravindran and van Neerven and therefore rule out the result of Altarelli and Lampe. 
However, we do have a theoretically-important difference compared to Ravindran 
and van Neerven, in that we find that the forward-backward asymmetry contains 
terms enhanced by logarithms of the quark mass. Even though these terms are 
numerically tiny for realistic quark masses, as a point of principle it means that the 
forward-backward asymmetry of massless quarks is not perturbatively calculable 
and non-perturbative fragmentation functions have to be introduced. 

We also calculated for the first time the corrections using the thrust axis 
definition rather than the quark direction. These lie approximately midway between 
the results of Refs. 0) and for the quark axis definition. 

We here only briefiy sketch the method and give the final result, and refer 
the reader to Ref. ® for more details. 



The simplest definition of tlie 6-quarkQ forward-backward asymmetry ApB 

is 

Nf - Nfi 

Afb = — — — , (1) 

wliere Np and Nb are tlie number of b quarks observed in the forward and backward 
hemispheres, respectively. 

The axis that identifies the forward direction can be defined in a variety 
of ways. In this paper we explicitly consider two different definitions: the 6-quark 
direction, and the thrust axis direction, which we denote by Ap^ and A^^^ respec- 
tively. 

According to the definition in Eq. (P, ApB can be expressed in an equiv- 
alent way in terms of the cross section 

da(e+e- b + X) 
dx dcos 9 

for inclusive 6-quark production, where x is the fraction of the electron energy carried 
by the b quark and 6 is the angle between the electron momentum and the direction 
defining the forward hemisphere (both energies and angles are defined in the centre- 
of-mass frame). 

Starting from the distribution in Eq. (|^), we can introduce the forward and 
backward cross sections ap and cr^: 

/■I dcr fO /■! dcr 

ap = / dcos 6* / dx- — -, aB = / dcos 6* / dx - — -, (3) 

Jo Jo dxdcost^ J-i Jo dxdcosty 

and the symmetric and antisymmetric cross sections as and a^'- 

o's = (^F + o-B, aA = CTp - aB ■ (4) 

We can then write the forward-backward asjTiimetry as 

ApB=^-^ . (5) 

In order to calculate this ratio perturbatively, we first separate the con- 
tributions to the cross sections into three classes: flavour non-singlet {NS), flavour 
singlet (S), and interference (or triangle) (Tr) (see Ref. i) for their precise defini- 
tion). We thus write the cross sections as 

as = (rs,NS + 4^s + 4!Tr + 0{al) , (6) 

(TA = (JA,NS + Cr^Tr + C(«s) • (7) 



^Throughout this paper we exphcitly consider the case of the 6-quark. The results for the charm 
quark can be simply obtained by properly replacing the mass and couplings of the massive quark. 



In this notation, up to 0{as) there are only non-singlet contributions. Thus, 
o's'^si '^^sj'r "^A.Tr proportioual to ctg. There are no singlet contributions 
to the antisymmetric cross section a^- 

The forward-backward asymmetry is decomposed in a similar way. Ex- 
panding the ratio a a/ 0-3 up to C(ag), we write 

(0) / (2) (2) (2) 

.(2) _ .(2) / ^A,Tr _ ^ S,Tr _ ^ S,S 



FB - ^FB,NS -r (q) (o) (0) (0) 

\ 



(2) 

where Ap^j^g denotes the non-singlet component: 

.(2) _ (^A,NS , . 

-^FB^NS — ■ V'l 

CS,NS 

The triangle contributions give non-universal (i.e. non-factorizable) correc- 
tions to both the symmetric and antisymmetric cross sections. They are calculated 
in Ref. 0) for the 6-quark axis definition and found to be very small. To our knowl- 
edge their contribution to the thrust axis definition has never been calculated, but 

(2) (2) 

we expect it to be similarly small. We therefore neglected it, i.e. (JsTr cr\ in 
Eq. (H), from our calculation. 

The singlet contribution to the symmetric cross section, as, is logarith- 
mically enhanced in the small-mass limit and proportional to agln^Q^/m^. An 
approximate expression for it, denoted by i^Branco^ ^^^g ^gg^ j^gf_ |^)_ jg calcu- 
lated exactly to 0(ag) in Refs. i' , and the leading and next-to- leading logarithms 
are summed to all orders in as in Ref. . 

In some sense the singlet component is a 'background' to the forward- 
backward asymmetry measurement and, in fact, in the experimental analyses (see 
e.g. Ref. 0)) it is statistically subtracted using Monte Carlo event generators. We 
therefore neglected it, i.e. in Eq. (H), from our calculation. 

' (2) 

Before describing the calculation of Ap^j^g, we take a slight diversion to 
discuss the contribution to it from four-6 final states. Let us first point out a basic 
fact. The four-6 process contributes to both the 6-quark cross sections as and a a 
and the e~^e~ total cross section. However, they appear with different multiplicity 
factors in the two cases. In the case of the e^e~ total cross section the multiplicity 
factor is simply equal to unity. In the contribution to the inclusive 6-quark cross 
sections as and cr^, these terms count twice since there are two b quarks in the final 
state. This observation is important in understanding the results for the non-singlet 
component of the symmetric cross section as discussed shortly. 

After summing and squaring the Feynman diagrams for four-6 production, 
we obtain two types of contribution: i) those that are identical to the bbqq final state 



but with the other quark q replaced by an untriggered-on b quark, and ii) those 
that are genuine interference terms arising from the fact that the two antiquarks 
are indistinguishable, called the ii^-term in Ref. 1^). The squared diagrams of type i) 
are lumped together with the corresponding terms from bbqq in the singlet 
in Eq. (|^)), non-singlet {as^NS and aA,NS in Eqs. and (J^)) or triangle {cr^sTr 
and o'A j'r in Eqs. and (0)) contributions. The squared diagrams of type ii), 
which give a universal (i.e. factorizable) correction to both the antisymmetric and 
symmetric cross sections, can be considered part of the non-singlet contributions. 

It is not entirely clear how four-quark final states are actually treated 
in the different experimental analyses, i.e. the extent to which they are genuinely 
measuring the inclusive cross sections. Often some vague statement like "a four-6 
final state is more likely to be tagged than a two-6 one, but less than twice as likely" 
is made. To know what to calculate one must understand the corrections that 
are applied for this difference in tagging efficiency, which are not usually explicitly 
stated in the papers. In the absence of a unique experimental procedure and of a 
definitive statement from the experiments on what they are measuring, we make this 
ambiguity explicit by multiplying the £^-term by an arbitrary weight factor We- An 
inclusive definition would correspond to We = 2 (each b quark contributing once), 
while an exclusive definition (the cross section for events containing at least one b 
quark) would correspond to We = 1. Since the forward-backward asymmetry is 
defined to be the asymmetry of a differential cross section, it is clear that we must 
use the same cross section definition in the numerator and denominator, i.e. that 
We must be the same in the symmetric and antisymmetric cross sections. 

Having defined the weight factor We for the ii^-term, we can define the 
following symmetric and antisymmetric cross sections 

as,NsiWE) = as,Ns{WE = 0) + We J Es , (10) 

<yA,Ns{WE) = (TaMWe = 0) + We J Ea , (11) 

where / Es and / Ea denote the integral of the symmetric and antisymmetric E- 
term, respectively. We recall that the 'truly' inclusive cross sections in Eq. (^ 
correspond to the definition with We = 2, i.e. crs,NS = crs,Ns(WE = 2) and <7a,ns = 

(TA,NsiWE = 2). 

The C(Q;g)-calculation of the cross sections in Eqs. (p!0| , |Tl|) and of the 
corresponding forward-backward asymmetry in the case of a finite 6-quark mass is 
extremely complicated, and we are not able to perform it. It is thus convenient to 
separate the calculation into a piece that is finite (although still cumbersome) in the 



massless limit and a simpler piece that is not. Then, the finite piece can be more 
easily computed in the massless approximation, while the simpler non-finite piece 
can be computed in the massive theory. 

It is possible to show @) that the inclusive definition, with We = 2, results 
in an antisymmetric cross section cr^ (or, analogously, crA,Ns) that is finite in the 
massless limit, at least at 0{al). However, in the same limit, the inclusive symmet- 
ric cross section is divergent at 0(«g), even if we only consider its non-singlet com- 
ponent. The corrections to (the non-singlet component of) the forward-backward 
asymmetry itself must therefore also be divergent in the massless limit. 

This final statement remains true for any value of We > 0. For example, 
with We = 1, the non-singlet part of the symmetric cross section is finite, but the 
antisymmetric cross section contains logarithmically-enhanced terms. 

The divergences in the non-singlet components correspond to logarithmi- 
cally-enhanced terms aglnQ^/mf coming from the ii^-term in the triple-collinear 
limit, i.e. when three fermions of the four-quark final state become simultaneously 
parallel. The integral of the symmetric E-teim is calculated numerically in Ref. @) 
and, neglecting corrections of 0{mb/Q), the final result is 



2 ( — - — + 2C3 ) In - 8.1790 ± 0.0013 



4 2 ^7 mi 

(12) 

where the analytic coefficient in front of In Q'^/m^ is proportional to the integral of 
the non-singlet Altarelli-Parisi probability P^^{z,as) (see, for instance, Ref. 1^)): 



/' 

Jo 



|)^c.(c.-lc.)(H_|^ 



dzP,f(^,as) = (;2) CE{CE--CA]{-i-- + 2C,\ . (13) 



and the constant term is the result of our numerical calculation. 

Having pointed out that the symmetric ii^-term is divergent in the massless 
limit, it is very simple to show how the divergence appears in the inclusive symmet- 
ric cross section. According to the definition of the non-singlet component of (T5, 
the virtual diagrams that contribute to cts^ns are exactly those that contribute to 
the e^e~ total cross section. As for the real diagrams, they only differ by the con- 
tributions coming from the E-ieim. In the total cross section, the E-ieim. enters 
with a multiplicity factor We = 1, and its divergence is cancelled by that of the 
virtual diagrams. In the inclusive 6-quark cross section, the multiplicity factor of 
the E-teim. is We = 2 and, thus, the cancellation of the divergence with the virtual 
terms is spoiled. 



This argument also allows us to directly compute the C(ag)-correction to 
Eq. (jlOl). Exploiting the fact that the massless QCD correction to <Js,ns(}Ve = 1) 
is equal to the correction -Rg+e- to the total cross section, we write 



c^s,NsiWE) = a 



(0) 



Re+e- + {We -1)JEs + 0{a, 



(14) 



Then, we obtain an explicit expression for (Js,NsiWE) by simply introducing in 
Eq. (|n|) our result in Eq. (|T^ for / Eg and the well-known result 0) for Re+e-- In 
particular, for the inclusive symmetric cross section we obtain 



(^s,NS = (^s,ns(}Ve = 2) = cr 



(0) 

s 



Re+e- + j Es + 0{al) 



(15) 



Since both aA,Ns{WE = 2) and as,Ns{WE = 1) are finite when 0, we 

can use the dependence on We to construct an unphysical observable that is finite 
in the massless limit: 

.(2);finite _ (^A,Ns{We = 2) 

= 77J7 TT ■ UOJ 

(yS,NS\WE = ij 

This observable is the ratio of the antisymmetric part of the inclusive cross section 
and the symmetric part of the exclusive cross section. Thus Ap^J^^^^'^ is unphysical 
in the sense that it is not the forward-backward asymmetry of a single differential 



cross section. Nonetheless the definition in Eq. (16) helps us to perform a massless 
calculation. The physical result for We = 2 is then given by 



.(2) 

^FB,NS 



^g),finite _ J (17) 



where / Eg is the integral of the symmetric i?-term, given in Eq. (12). 



Even in the massless limit numerical two-loop calculations are prohibitively 
difficult to set up. Fortunately there is a cancellation between the genuinely two- 
loop effects in the ratio on the right-hand-side of Eq. ([I6|), which does allow its 
numerical evaluation. The total contribution can be written as 

n ■ _L ^(l);one-loop (l);tree (2);two-loop (2);one-loop (2);tree/... _ 

^(2);finite _ a + + A + cr A + cr^ [We - ^) 

FB - ^^0) ^ ^a);one-loop ^ ^a);tree ^ ^^2);two-loop ^ ^^2);one-loop ^ ^imree^^^ ^ 

(18) 

The 0(Q;s)-contributions come from the one-loop cross sections a^i^'O'i^-looP for the 
two-parton process e"'"e~ bb and the tree-level cross sections a^^"^'^^^^ for the three- 
parton process e+e~ — > bbg. Similarly the non-singlet (9(ag)-contributions from the 
two-parton, three-parton and four-parton final states are denoted by (j(2);two-loop^ 
^(2);one-loop g^-^^ ^(2);tree respectively. Of course, the dependence on We enters only 
through the four-parton terms a''2^'^^'^^{WE = 2) and af^'^^'^'^{WE = !)• 



Each of the cross sections is separately divergent, so they have to be regu- 
larized in some way before being combined together. In any regularization scheme 
that preserves the helicity conservation of massless QCD0 (for example, dimensional 
regularization), we have the properties 



(l);one-loop (l);one-loop (2);two-loop (2);two-loop 

^ ^^4 ^ 

^(0) _(0) ' (0) (0) 



(19) 



so that if we expand the ratio in Eq. (|T^) up to 0(ag), the two-loop corrections 
cancel, and we obtain 



(0) 

.(2);finite _ (^A 
^FB — ~ 



a 



's 

a 



(0) 




O c 



^(0) / I (0) (0) 
^ ^^A ■ 




(20) 



(2);one-loop (2);one-loop (2);tree/Tj^ on (2);tree/Tjr n 
a\" ^ a);'' a\" {We = 2) ay' [We = I) 



JO) JO) ' (0) (0) 

cr^ (T^ ag 

where a^^^ and cr^ ■* are the complete contributions to the antisymmetric and sym- 
metric cross sections at 0{as), o"^^-* = (j(i);one-loop _j_ ^{i);tree_ rpj^g gj,g^ jj^^g ^^-^ 

calculated analytically, but the second line is too comphcated to be able to, so must 
be done numerically. Since the two-loop terms have cancelled, this has the structure 
of a NLO three-jet calculation, as first noticed by Altarelli and Lampel^). Thus the 
calculation can be performed using known techniques (we use the dipole-formalism 
version of the subtraction method 

We are finally ready to present our numerical results. We start with the 
unphysical, but finite, quantity defined in Eq. (|l^), and separate out the different 
colour factors, as in Refs. i D: 



-(2);finite;b _ .(0) 
^FB — ^FB 



^^-Ce) (-Ce 
271 V 27r2 / V2 

+ (^) ' Ce {CCe + NNc + TTnNf] 



(21) 



with as = <ys{Q^)- Our numerical results are shown in Table |l|, in comparison 
with the previous calculations. It is clear that we disagree badly with the results 
of Altarelli and Lampe ^ , but are in excellent agreement with Ravindran and van 
Neerven 0), who give the coefficients analytically. However, we should recall that 
Ap^J^^^^'^, as given in Eq. (pi]), is not the forward-backward asymmetry of a definite 
cross section. The physical forward-backward asymmetry must have subtracted 
^Note that the relations (p^) are expHcitly violated for massive quarks. 



6-quark axis 


C 




T 


AL a 
RvN D 


4.4 ±0.5 
\ = 0.375 


-10.3 ±0.3 
_1|3 ^ _i5.375 


5.68 ±0.04 
^ = 5.5 


Our Calculation 


0.3765 ± 0.0038 


-15.3769 ±0.0034 


5.5002 ±0.0008 



Table 1: Results for the coefficients of the 0(as) correction to the finite part of the 
forward-backward asymmetry with the b-quark axis definition, Eqs. ( fr?| , W^) - 



from Eq. (|2T|) the logarithmically-enhanced term of Eq. (17), which is not present in 
the result of Ref . ® . Thus it seems that they have somehow computed the unphysical 
^FB^'^^*^ rather than the forward-backward asymmetry. In fact, their expression for 
the correction to the symmetric cross section (/t±/l in their Eqs. (31) and (32)) is 
actually equal to our (Js,ns{We = 1) in Eq. |l^. So, the fact that their result for A^p\ 
agrees with our A^^^'^^*^ means that we confirm their result 0) for the inclusive 
antisymmetric cross section cr^-* = (t^^{We = 2) (/^ in Eq. (33) of Ref. i)). 

Using our numerical program it is straightforward to calculate the forward- 
backward asymmetry with any other axis definition (or cuts, for example on the value 
of the thrust). With the thrust axis definition, we obtain 



i (2);finite;T 



A 



(0) 
FB 



1 - 



OS 

27r 



± 



as 3 
' 27r2 

27r/ 



Cpj {IMCp) 

Cf {CCf + NNc + TTnNf) 



(22) 



with as = «s(Q^) and the coefficients given in Table 0. The logarithmically- 
enhanced piece that has to be added to this is identical to that in the 6-quark 
axis definition, namely Eqs. (|1^, |T^). It is worth noting that the difference between 
the two definitions is the same size and in the same direction as at 0{as), leading 
to an overall difference of 0.8% for as ~ 0.12. 

We finally recall that we include an arbitrary factor We in front of the 



thrust axis 


C 


N 


T 


Our Calculation 


-3.7212 ±0.0065 


-9.6011 ±0.0049 


4.4144 ±0.0006 



Table 2: Results for the coefficients of the 0(ag) correction to the finite part of the 
forward-backward asymmetry with the thrust axis definition, Eqs. ( |7^ , ^Bj ). 



four-6 contribution to account for the way in which it is treated in the experimental 
analyses. For a fully inclusive definition, in which each h quark contributes once, 
We should be set equal to 2, while for an exclusive definition. We should be set 
equal to 1. Our final result for the non-singlet component of the forward-backward 
asymmetry, is then: 

(yA,Ns{WE) 



(yS,Ns{WE) 



.(2);finite .(0) 



FB 



^We) ( 2 



Ea - 



Es ]+lWEj Es 
(23) 



where Ap^J^^^^*^ is given in Eqs. (|2T|, |2^) and Tables |l| and ^, / Es is given in Eq. 
and (see Appendix B of Ref . ® ) 



Ea- 
Ea- 



Es 
Es 



as 
2n 



Cf{Cf - iC^) (0.3620 ± 0.0007 

0009 



JEa-JEs = (g) C^(C^-iC^) (0.1144 ±0. 



quark axis, (24) 
thrust axis. (25) 



Note that the combinations of E'-term contributions in Eqs. 
in the massless limit (see the discussion in Appendix B of Ref. 

Putting all these numbers together, and setting Nj = 
forward-backward asymmetry according to the two definitions as: 



and (|25D are finite 
5, we write the 



FB,NS 



A 



(2);T 
FB,NS 



(We 



(We) 



A 



(0) 
FB 



0.318as - 0.973a^ + WEal fo.004051n^ - 0.0240 j 



(26) 



^(0) 
^FB 



1 - 0.284as - 0.676a| + WeoI |^0.00405 In - 0.0233 j 



(27) 



with as = cts(Q^)- Note that the logarithmically-enhanced term, InQ^/m^, is 
present for any physical {We > 0) value of We- 

Putting in an explicit value for as, we summarize the total QCD correction 
according to the various available calculations in Table ^. We continue to neglect all 
terms that vanish in the massless limit. Since in the existing experimental analyses 
(see for example Ref. 0)), the known 0{as) correction for the thrust axis definition 
was included, together with the Altarelli and Lampe quark axis value for the 0{ag) 
corrections, we do the same in Table |^. 

We find that the difference between the Ravindran and van Neerven calcu- 
lation and ours is numerically irrelevant, being smaller than 10~^ for b quarks and 
~ 2.5x10"^ for c quarks. Therefore at the numerical precision required by current 
or any foreseen experiments, we agree with their result - the difference is only one 





AL 
quark axis 


RvN 1) 
quark axis 


Our Calculation 
quark axis 


Our Calculation 
thrust axis 


Correction, Ap^^/Ap^^ 


0.962 


0.952 


0.952 


0.956 



Table 3: Total QCD correction to the forward-backward asymmetry in the small- 
mass limit, with as = 0.12. In each case, the thrust axis definition is used for 
the 0{as) correction and the definition shown is used for the 0{a'g) correction, as 
discussed in the text. 



of principle. The difference between the Altarelli and Lampe calculation and ours 
for the quark axis definition is more significant though, at around 1%. However, the 
error in their calculation and the effect of using the thrust axis definition partially 
cancel, and the total difference is around 0.6%. 

We should also mention the important fact, discussed in Ref. 0), that the 
experimental procedures introduce a bias towards more two-jet-like events. This 
actually decreases the size of the QCD corrections considerably, so our number 
should be considered as an upper bound on the final difference. 

In Ref. @) we try to estimate the remaining uncertainties in the forward- 
backward asymmetry, bearing in mind that while the 2% precision of current exper- 
iments is close to their final limit, a future linear collider could be capable of experi- 
mental errors on the left-right forward-backward asymmetry of order 0.1% HD. We 
found several sources of uncertainty that all contribute at the few per mille level. 
While this is certainly sufficient for the current precision of the data, matching the 
precision of a future linear collider measurement could be extremely difficult. It is 
likely that this could only be done by making even more stringent two-jet cuts in 
order to work in a region in which the corrections and their uncertainties are smaller. 
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